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ABSTRACT 

In this model we study the dynamics and control of 
hepatitis B virus (HBV) infection which is a major 
health problem worldwide by considering condom, 
vaccination and treatment as control measures. 
Initially we determined the basic reproduction number 
R_0 for the model and observe that once R_0<1, the 
disease free equilibrium will be stable and HBV 
infection can be controlled using the three control 
measures and we also study the solution of the 
endemic equilibrium point of the model. Next we take 
the sensitivity analysis of the basic reproduction 
number of HBV infection and obtain that the 
endemicity of the infection will reduce with the 
controls. Finally, the numerical simulation result 
shows that combination of condom, vaccination and 
treatment is the most effective way to control hepatitis 
B infection. 

Keywords: Mathematical model, HBV Infection, HBV 
Control, Condom, Vaccination, Treatment. 

I. INTRODUCTION 

Hepatitis B was called serum hepatitis in the past 
(Southern Nevada Health District (SNHD), 2015). 
This virus has a number of serological types and the 
infection occurs at any age but mostly in adults. The 
incubation period is 50 to 180 days. The virus enters 
the blood and is spread by blood and blood products. 
Antibodies are formed and immunity persists after 
recovery. Infection usually leads to severe illness 
lasting 2 to 3 weeks, often followed by a protracted 
convalescence. Carriers may or may not have had 
clinical disease. Type B virus (serum hepatitis) may 
cause massive liver necrosis and death and in less 
severe cases recovery may be complete (Corazon, 


2015). In chronic hepatitis which may develop, live 
virus continue to circulate in the blood and other body 
fluids, that is, saliva, semen, vaginal secretions. The 
condition may predispose to liver cancer. Hepatitis is 
the inflammation of the liver caused by virus, toxic 
substances or immunological abnormalities. 
Infectious hepatitis is caused by viruses, several types 
of which have been isolated as specific causes of the 
disease and can be detected by tests. 

Associated with acute viral hepatitis is acute infection 
with hepatitis B virus and it starts with general ill - 
health such as mild fever, dark urine, nausea, loss of 
appetite, vomiting, body aches, and then progress to 
development of jaundice. The illness lasts for few 
weeks before gradually improves in most effected 
people. Few patients may die as a result of fulminant 
hepatic failure and in some individuals, the infection 
may go unrecognized (asymptomatic) (Pan and Zeng, 
2005). Chronic HBV may be either asymptomatic or 
associated with a chronic inflammation of the liver, 
leading to cirrhosis over a period of several years. 
Incidence of hepatocellular carcinoma (liver cancer) 
increases in chronic HBV patients. Avoidance of 
alcohol consumption is recommended to chronic 
carriers of HBV infection because it increases their 
risk for cirrhosis and liver cancer. Hepatitis B virus 
has been linked to the development of membranous 
glomerulonephritis (MGN) (Chang, 2007). 

II. Related Literature 

Hepatitis B virus infection can be controlled either 
through prevention or treatment. The most effective 
way of preventing HBV infection is by immunization 
with hepatitis B vaccine. Prevention of HBV 
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transmission during early childhood is important 
because of the high likelihood of chronic infection 
and chronic liver diseases that occur when children 
less than five years of age become infected. In HBV 
infection, the development of anti - HBs as HBsAg 
disappears suggests that anti - HBs are the protective 
antibody. Therefore, immunoprophylaxis of HBV 
requires the administration of anti - HBs antibodies 
by vaccination with HBsAg. The vaccine schedule 
most often used for adult and children has been three 
dose intra muscular injections in which the second 
and third doses are administered 1 and 6 months after 
the first. Recombivax HB has been approved for aged 
11-15 years. In 2007, Twinrix was approved as a 
four dose accelerated schedule by Food and Drug 
Administration (FDA) and the schedule of doses are 
0, 7, and 21 - 31 days and a booster dose at 12 
months after the initial doses. If the interval between 
doses is longer than the recommended interval, it is 
not necessary to restart the series or add doses (Center 
for Disease Control and Prevention (CDC), 2012). 
This vaccine is for persons at the age of 18 years and 
above. Also, there is vaccine for infants born to 
HBsAg negative mothers (coffin et al, 2011). The first 
dose can be administered before the newborn is 
discharged from the hospital. The second dose can be 
given after 1 to 2 months and the third dose can be 
given 6 to 18 months later. One cannot be infected 
with hepatitis B from these vaccines (Chu and Liaw, 
2007). Following vaccination, hepatitis B surface 
antigen may be detected in serum for several days; 
this is known as vaccine antigenaemia (Arauz-Ruiz et 
al, 2002). Protection has be observed to last 12 years 
in individuals who show adequate initial response to 
the primary course of vaccinations, and that immunity 
is predicted to last at least 25 years (Glebe and Urban, 
2007). 

HBV infection can also be prevented through the use 
of condom or another type of latex barrier while 
having sex because it have been shown that consistent 
and correct use of condom is 98 - 99% effective in 
preventing HBV transmission. Since HBV infection is 
transmitted through body fluids; prevention is thus the 
avoidance of such transmission like unprotected 
sexual contact (use of condom), blood transfusion, 
reuse of contaminated needles and syringes, and 
vertical transmission during birth. Acute hepatitis B 
infection does not usually require treatment because 
about 95% of most adults clear the infection 
spontaneously (Shapiro, 1993, Lok et al, 2015). Early 
antiviral treatment may only be required in fewer than 


one percent of patients, whose infection are of great 
severity (fulminant hepatitis) or who are 
immunocompromised. In patients with fulminant 
infection, intravenous (IV) feeding may be needed, if 
the patient has persistent vomiting and can’t take oral 
foods. Ninety percent of patients with acute HBV 
infection have a favorable course and recover 
completely (Jayapal, 2007). On the other hand, 
treatment of chronic HBV infection may be necessary 
to reduce the risk of cirrhosis and liver cancer. 
Chronically infected individuals with persistently 
elevated serum alamine aminotranferase, a marker of 
liver damage and HBV DNA level are candidates for 
therapy (Shi et al, 2011). Thirty to forty percent of 
patients with chronic hepatitis B respond to IFNa 
treatment. 

Immune defense is partially active in chronic hepatitis 
B patients, since the viral load and HBsAg level are 
lower (Kosinska et al 2013). Undertaken in various 
experimental systems is antigen expressing DNA in 
combination with antiviral therapy or immunotherapy 
with various HBV derived antigens and example is 
the work of Michael Roggendorf in infected 
woodchuck with partial success (Kosinska et al 2013). 
Hepadnaviral cccNDA silencing or destruction would 
be the best therapy for control of HBV infection but it 
is still a science fiction (Gerlich, 2013). The current 
therapy (GS - 9620) for HBV infection suppresses 
viral replication and delay progression of liver 
damage. This GS - 9620 which works by targeting the 
virus represents the first conceptually new treatment 
for hepatitis B virus in more than a decade (Texas 
Biomedical Research Institute, 2013). Other antiviral 
medicine includes, Lemivudine (Epivir) and the 
newer neucleotide analogues famciclovir, lobucavir 
and adefovir (Hepsera), dipivoxil, tenofovir (viread), 
tesbivudine (Tyzeka), entecavir (Baraclute) (Lok et al, 
2015) and they directly block or suppress or destroy 
the HBV by interfering with its replication (Gerlich, 
2013; WHO, 2015). These agents are highly effective 
and well tolerated. IFNa and lamivudine have been 
approved by FDA for the treating chronic hepatitis B. 
For decompensated liver disease, liver transplant may 
be required, i.e., surgery may be required to treat 
hepatocellular carcinoma (Bonino, 1987and Jayapal, 
2007). However some individuals are much more 
likely to respond than others and this might be 
because of the genotype of the infection virus or the 
patient’s heredity. The treatment reduces viral 
replication in the liver, thereby reducing the viral load 
(the amount of virus pareticles as measured in the 
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blood)(Lok and Mcmahon, 2007). Response to 
treatment differs between genotypes. Interferon 
treatment may produce an e antigen seroconversion 
rate of 37% in genotype A but only a 6% 
seroconversion rate in type D. Genotype B has similar 
seroconversion rates to type A (Zoulim, 2004). 

Garib - Zahidi (2012), used the model presented by 
Nowak et al (1996) to study the kinetics of HBV 
infection (cellular automation). They showed that 
cellular automation is sophisticated enough to 
reproduce the basic dynamical features of the 
infection. They acknowledged that the basic 
dynamical features of the infection are the biphasic 
behavior of viral load and age of various patients. 
Kamyad et al, (2014) constructed a model with two 
controls (vaccination and treatment) based on the 
characteristics of HBV transmission. They observed 
from the numerical results that mixed control 
measures responds better than any other independent 
control. 

Abdulrahman et al, (2013) in a model of HBV 
transmission dynamics in a population with viral 
dynamics, incorporating vertical transmission and 
sexual maturity showed that all the parameters in their 
model are sensitive to the transmission and prevalence 
of HBV either positively or negatively and that the 
most influential of the parameters are the natural 
death rate (pi), birth rate (b), HBV - sexual 
transmission probability contact rate (p), and the total 
sexual contact rate (c), with the last two parameters 
having +1 sensitivity index.In the analysis of the 
transmission dynamics and control of hepatitis B virus 
(HBV) infection in Nigeria using six different control 
strategies with a 25 years waning rate of vaccine, it 
was observed that the present stratege of HBV 
infection control can never bring about disease control 
in Nigeria (Abdulrahman et al, 2013). They also 
stated that effort should be made so that 20% to 30% 
of sexually active susceptible Nigerians will be 
vaccinated yearly, to be added to the 66% who are 
vaccinated yearly at birth. 

Momoh et al (2012) confirmed that much progress 
have been made in understanding the transmission 
dynamics of HBV infection but maintained that more 
effort is needed if we want to conquer the infection. 
Obaid and Elaiw (2014) proposed and analyzed two 
virus infection models with antibody immune 
response and chronically infected cells. It is not 
reachable in practice to vaccinate newborns by a rate 


greater than 96%. Therefore, treatment or additional 
vaccination strategies are needed to control the spread 
of HBV infection in population (Moneim and Khalil, 
2015). Certainly, controlling HBV infection simply 
indicates reducing the susceptible by vaccination or 
reducing the latent and infected by treatment scheme 
Kimbir et al (2014), extended the model of Zou et al 
(2009) and used it to study the transmission dynamics 
of HBV infection considering vaccination and 
treatment as control measure in the host population. 
They showed from the numerical simulation of their 
model that effective vaccination and treatment is a 
good control strategy for HBV infection. But a 
combination of vaccination and treatment is a better 
control strategy for the disease. In the recovery of 
HBV infection, each component of innate and 
adaptive immune response contributes, but in the 
absence of one component of innate immunity, the 
remaining two defense mechanisms are sufficient for 
viral clearance (Dontwi et al, 2010). 

III. MODEL FORMULATION 

3.1: THE ASSUMPTIONS OF THE MODEL: 

1. We assume that both Susceptible, Exposed 
(infected), infectious and Recovered individuals 
die at the same natural rate pi. 

2. The treated individuals recover, unless the 
individual dies naturally (Ahmed et al, 1987 and 
Oleary et al, 2008). 

3. Also, we assume that vaccinated individuals do 
not acquire permanent immunity because the 
efficacy of the vaccine wanes (Keeling and 
Rohani 2008). 

4. Exposed (infected) individual can spontaneously 
recover from the infection 

5. Infectious but not treated individuals can recover 
by the effect of the natural immune response 
(noncytocidal responses) (Zou et al, 2010; Kar and 
Jana 2013). 

6. Once an individual progresses to infectious class, 
he is classified as infectious but treated or 
infectious but not treated. 

3.2: MODEL VARIABLES FOR CONTROL 

OF HEPATITIS B VIRUS INFECTION 

The equations are partitioned into twelve 
compartments of: 

S(t) Susceptible individual; 

S v (t) Susceptible but vaccinated; 

S u (t) Susceptible but not vaccinated; 

S vc (t) Susceptible, vaccinated and uses condom; 
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S uc (0 Susceptible, not vaccinated but uses 

condom; 

S V nc (0 Susceptible, vaccinated but does not use 
condom; 

Sunc (0 Susceptible, not vaccinated and does not 

use condom; 

E(t ) Exposed (infected); 

E r (t) Exposed (infected) and treated; 

7 r (t) Infectious and treated; 

/ nT (t) Infectious and not treated; 

R(t) Recovered. 

3.3: MODEL FORMULATION FOR HBV 

CONTROL 

The model for control of HBV in human population 
were formulated using xp as the recruitment level of 
the population at time t. This model incorporated 
vaccination given to susceptible individuals, condom 
use by vaccinated and non vaccinated individuals and 
treatment given to exposed (infected) and infectious 
individuals. In the model, we made use of the 
following variables and parameters. 

3.4: MODEL PARAMETER FOR HBV 

CONTROL 

xp 1 Recruitment level of susceptible 

individuals; 

p 1 Rate at which susceptible individuals are 

vaccinated; 

a(l — p x ) Rate at which susceptible individuals are 
not vaccinated; 

/?! Rate at which vaccinated susceptible 

individuals use condom; 

a 5 ( 1 —/? x ) Rate at which vaccinated susceptible 
individuals does not use condom; 
a Denote the rate of condom usage by non- 

vaccinated individuals; 

b(l — a) Denote the rate at which susceptible 
individuals who are not vaccinated do 
not use condom; 

0 1 t 1 The rate at which susceptible individuals 

who are vaccinated but does not use 
condom are being exposed to the 
infection due to waning of the vaccine; 
5(1 — c) The rate at which non-vaccinated 
susceptible individuals who use condom 
are being exposed to the infection due to 
error in condom/condom usage; 
v Denote the rate constant at which 

susceptible individuals who are not 
vaccinated and does not use condom 
move to exposed class 


a i 


«2 


1 — a 1 


/c 4 ( 1 - d) 


m 


Pi 


y i 


9i 


Pi 


The rate of exposed (infected) individuals 
who obtain treatment; 

The rate of exposed (infected) individuals 
who had spontaneous recovery; 

The rate at which exposed individuals do 
not obtain treatment and moves to 
infectious class; 

The rate at which infectious individuals 
are treated; 

The rate at which infectious individuals 
are not treated; 

The rate at which treated infectious 
individuals move to recovered class; 

Rate at which non treated infectious 
individuals move to recovered class; 
Death rate of non treated infectious 
individuals due to the infection; 

The rate at which treated exposed 
individuals move to susceptible class; 

The rate at which recovered individuals 
move to susceptible class; 

The movement rate of vaccinated 
susceptible individuals who use condom, 
to the susceptible class; 

The natural death rate of both the 
susceptible class, exposed (infected) 
class, infectious class and Recovered 
class. 


Wi 


Yi 



Fig 1: Flow diagram of Hepatitis B Virus Control 
Model 
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From the diagram above we obtain the following equations; 
dS 

= V'i - Pi s ~ a (! - Pi)S ~ Pi s + e i s vc + Pi^Y + YiR 
dS v 

— = PiS - Pis v - /?! s v - a 5 (l - p i )S v 


dt 


dS, 


dt 


r = a( 1 - Pi)5 - PtS u -aS u -b( 1 - cr)S u 




dt 


P\Sy Pl^VC 




dt 


Pl^UC 5(1 


dS n 


dt 
dS ,, 


&5(1 Pi)S v Pi^vnc 0lTl*-*u 
— 5(1 Pl^unc ~ vS U nc 


(3.1) 


dE 

dt 


dt 

tpiT-i^vnc T 5(1 c)S uc + vS unc (1 oc^E (XiE PiE 

dE r 


dt 


— a 1 E — PxEj — PiE T 


dl 


nT 


dt 


dI T 

-j-j- — d(l — a^I-p — — ^I-p 

— ^ 4(1 — 5)(1 — Cti)/ n r — tn/ nr — PlInT — ylnT 


dR 

— = ^I T + yl nT - 1 - a 2 E - p x R - YiR 


3.5: ANALYSIS OF THE SYSTEM FOR free equilibrium point and the endemic equilibrium 

CONTROL point 

If we assume that the control parameters 

Pi(t),^ 1 (t),cr(t),a 1 (t),and d(t) are constant 3.7: DISEASE FREE EQUILIBRIUM 
functions, we analyze the equilibrium point and the The disease free equilibrium state denotes absence of 
basic reproduction number of the system. infection, vaccination, use of condom and treatment. 

Therefore, all the infected class will be zero and the 
3.6: EQUILIBRIUM POINT AND BASIC entire population will comprise of Susceptible, 


REPRODUCTION NUMBER 


Vaccinated, Non vaccinated, Condom compliance and 


We shall discuss the two possible non negative Non condom compliance individuals. At equilibrium 
equilibrium of system (3.1) and these are the disease state, we set the rate of change of each variable to be 

equal to zero. 

Thus 


dt dt 


dt 


dt 


dt 


dt 


dt 


_ I _ “'ll) _ ““ _ Q 

dt dt dt dt dt 


(3.2) 


AlsO let, (S, S v , S u , S vc , S uc , S vnc , S unc , E, E T , I T , I nT , R) — (S , S v , S u , S vc ,S' uc ,S vnc ,S unc , E ,E T , I nT , R ) at 
equilibrium state. 
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Therefore, 

0 = ipi — Pl S — a( 1 — Pl )S — Pl S + 0iS vc + Pl E T + YiR 
0 = PiS — Mi S v - /?i S v - a 5 (l - 
0 = a(l - Pl )S - Pl S u - oS u - b( 1 - a)S u 
0 Ml Syc dlSvC 

0 El^yc 5(1 ^)S UC 

0 ^ 5(1 Pi)S v Ml ^vnc ^Pl^l^vnc 

0 h(l 0")S U Ml Sync vSync 

0 = 01 TiSync + 5(1 - c)S uc + vS unc - (1 - aJE - a x E - a 2 E - Pl E 
0 = a t E — Mi E t — Pl E T 
0 = d(l — cc-l^It — Mi It ~ %ih 
0 = fc 4 (l - d)( 1 - ai)/ nT - m/ nr - MiW - T4r 
0 = ^/ r + y/ nT + a 2 F - MiR - YiR 


(3.3) 


Let fc s = a(l — Pl ), k 6 — a 5 (l — /?!), k 7 — 5(1 — c), k 8 — 5(1 — cr), k 9 — d( 1 — a!), 

k 10 — k 4 (l — d ). If we label the equations in (3.3) as (i) to (xii) and substitute E* — Ej — lj— I^ T — 0 into 
(3.3), we have from (3.3xii) that 

Zih + ylm ~ Ei r - YiR = 0 
(Mi + Yi)R = 0 
R = 0 


From (3.3iv), we have 


From (3.3ii) we obtain 


PiS v - P S VC - 9 ^ = 0 

PlSv 


s — 

J vc 


PlS 


Mi + 0i 


V Pi + Pi + ^6 


(3.6) 


(3.4) 


(3.5) 


Substituting (3.4), (3.5) and (3.6) into (3.3i) we get 

V>i - PiS - a(l - p,)S - Pi S + 0i (--4-) ( 0 PlS ) + 0 + 0 = 0 

VMi + ^i' \Pi + Pi + K-6' 

0i (Mi + 0i)(Mi + Pi + K) 


S = 


(Pi + Pi + k 5 )(Mi + 0i)(Mi + Pi + k 6 ) - 0i/?iPi 


(3.7) 


From (3.6) we have that 


5 - 

V Pl + /?1 + 


If we substitute (3.7) into (3.6) we obtain 
S„ = 


Pi0i(Pi + 0i) 


V (Pi + Pi + ^ 5 )(Pi + 0i)(Pi + Pi + k 6 ) - 6 i/?iPi 
From (3.3iii) we have that 

a(l - Pi)S - Pl S u -aS u -b(l- a)S u = 0 
(Mi + cr + k 8 )S u = k 5 S 
k 5 S 

S — --- 

u Pl + a + k 8 


(3.8) 


(3.9) 
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0i(Pi + 0i)Oi + Pi +k 6 ) \ 


Substituting (3.7) into (3.9) gives us 

^5 ( 


S„ = 


11 + a + k 8 V (Pi + Mi + ^ 5 ) (Mi + 0i)(Mi + Pi + k 6 ) — Q1P1P1, 


(3.10) 


Remember that (3.5) is 


S — 

^ VC 


Pi S v 


Pi + #1 


Therefore, if we substitute (3.8) in (3.5) we obtain 

S = ( _ h£i _)(_ 

VC HMi + #i)(Mi + Pi + k 6 )J \(Pi 


\p(n + di)(fi + pi + k 3 ) 


+ Mi + ^ 5 ) (Mi + #i)(Mi + Pi + — 9 1 (lp 1/ 

Pl^PlPl 

r P (3.11) 


VC (Pi + Mi + k 5 )(Mi + PiXMi + Pi + k 6 ) — QiPPi 


From (3.3v) we have that 


aS u - MiS wc - 5(1 - c)S uc = 0 


S — 


(tS-,, 


Mi +k 7 


(3.12) 


Then, substituting (3.10) into (3.12) we get 

ak 5 \f ^ 1(91 + 0i)(Pi+ Pi+ k 6 ) 


j) UPi 


UC ' (Mi + k 7 )(Mi + o + k 8 )J \ (p x + Mi + ^ 5 )(Mi + #i)(Mi + /?i + ^ 6 ) — 0i/?iPi, 


(3.13) 


Notice that (3.3vi) is 


S — 
°vnc 


£[ 5(1 PP)S V Pi^vnc 01^1 Sync 0 

(Ml 3” 01^1 Xvnc kfrS v 
k^S v 


Pi + 0 1 T 1 


(3.14) 


We then substitute (3.8) into (3.14) to obtain 

k 6 


S — 
•-’vnc 


s — 

•-’vnc 


Pi + 01 
^6pl 


-) ( _ el _ ) (— 

tJ Vmi + Pi + K> V(Pi 

■ J ((pi 


0i(Mi + 0i)(Mi + 0i + ^6) 


+ Mi + ^ 5 ) (Mi + #iXMi + Pi + k 6 ) 
0i (Mi + 0i) \ 


VMi + 01^7 V(Pi + Mi + k 5 )(pi + 0x)(mi + Pi + k 6 ) - fl^iPi, 
Also, we have that (3.3vii) is 


diPiPij 

(3.15) 


c — 

^unc 


^(1 Pl^unc V^unc 0 

(Ml 3” V^Sunc — k 8 S u 


Mi 3-u 

If (3.10) is substituted into (3.16), we obtain 
Mi + v Vmi 3- er + /c 8 ) ((p x 


(3.16) 


c — 

^unc 


0i (Mi 3- Si) (Mi 3-0i 3-^ 6 ) 


3- Mi + ^sXMi + 0i)(mi + P 1 + k 6 ) — 6 i p i p 1/ 


(3.17) 


Let Eq be the disease free equilibrium point, then, 

E* 0 = [S,S v> S u ,S vc> S uc> S vnc> S unc ,0,0,0,0] 
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Using (4.60) and substituting the value of k 6 and k 5 we have 

c _ 0i(Pi + 0i)(Pi +Pi + k 6 ) 

•JQ — 


•Xo — 


(Pi + Pi + ^5) (Pi + ^iXPl + Pi + k 6 ) — @lPlPl 
^i(Pi + fli)(Pi + q 5 ) _ 0i 

(Pi + a)(Pi + 0i)(Pi + a 5 ) - 0 p 1 + a 
_Pi0i(Pi + fli)_ 

(Pi + Pi + k 5 )(pi + Pi)(Pi + Pi + k 6 ) — 6 1 p 1 p 1 

(0)^i (Pi + Pi) 


(Pi + a ) (Pi + Pi)(Pi + a 5 ) ~ 0 


= 0 


$uo — 


0i(Pi + 0i)(Pi +/?i + k 6 ) 


Pi + o + fc 8 V(pr + Pi + /c 5 )(P + Pi)(P + /?i + h) ~ QiPiPi, 


a 


0i(Pi + Pi)(Pi + « 5 ) 


axp 1 


Pi + b \(Pi + a)(p x + PJGu + a 5 ) - 0/ (p x + a)^ + h) 

0i/?iPi 


SveO 


(Pi + Pi + k 5 )(pi + Pi)(Pi + /?i + k 6 ) - 
0i(O)(O) 


(Pi + a)(p! + PJCpi + a 2 ) - 0 


- 0 


c — 
‘-’ucO — 


akr 


0i(Pi + Pi)(Pi + /?i + k 6 ) 


(Pi + ky)(Pi + cr + P 8 )/ \(Pi + Pi + k s )(P + Pi)(P + /?i + k 6 ) - 9 1 f] 1 p 1/ 


(0)a 


0i(Pi + Pi)(Pi + a 5 ) 


(pi + P 4 )(Pi + b)J \ (pi + a)(p! + PJCp! + a 5 ) - 0 . 


= 0 


} vncO ~ 


k sPl 


Pi + 01 

' a 5 (0) 


u) ((Pi 


0i (Pi + Pi) 


+ Pi + /c 2 )(P + Pi)(p + /?i + ^ 3 ) — Pl/^lPly 
0 (p + Pi) 


,pi + 4> 1 t 1 \(pi + a)(p x + PJCp! + a 5 ) - 0 , 


- 0 


J uncO 


( —-— U — 
Pi + u Vp! + a + /c 8 / VCpi 


0i (Pi + Pi)(Pi + /?i + ^ 6 ) 


+ Pi + k 5 )(p + Pi)(p + /?i + k 6 ) — Pi/?iPi 


a 


(Pi + Pi)(Pi + « 5 ) 


abxp 1 


Pi Vp! + b/ \(Pi + a) (p x + PJCpi + a 5 ) - 0y Pi(pi + a)(p x + b) 

The disease free equilibrium point of the model Eq = [S’, S v , S u , S vc , S uc , S vnc , S unc , 0,0,0,0] is 

dh- 0 _^_ 0 0 0 _ ^ _0 0 ool 

7 H+a’ ’ (/ii+a)(/ii+a)' ' ' ViOi+OOi+a) ’ ' ' ' T 


3.8: STABILITY OF THE DISEASE FREE EQUILIBRIUM (DFE) 

In equation (3.1), all the parameters are non negative and we assume stable population with per capita 
recruitment of susceptible individuals and death rate (both natural and disease induced). For simplicity, the 
population size is normalized to one, that is S, S v , S u , S vc , S uc , 

S V nc> S unc , E, E r , I T , I nT , R are respectively the fraction of the susceptible, susceptible with vaccination, 
susceptible without vaccination , susceptible with vaccination and uses condom, susceptible with vaccination 
but does not use condom, susceptible without vaccination but uses condom, susceptible without vaccination and 
does not use condom, exposed, exposed but treated, infectious but treated, infectious but not treated and 
recovered. This can be represented thus; 

S + S v + S u + S vc + S uc -I- S vnc + S unc + E + E T + I T + I nT + R — 1 (3.18) 


According to Pang, (2010) and Kamyard, (2014) equation (3.2) holds and we can omit the equation for R in our 
analysis so that (3.1) becomes 
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dS 

fa = 0i ~ Pi s ~ a(l ~ pi)S ~ PiS + diSvc + PiE t 
+Ti[1 — (S + S v + S u + S vc + S uc -I- S vnc -I- S unc + E + E t + I T + I nT )] 




dt 


— PiS p,S v Pi S v <^ 5(1 Pi)S v 


dt 


= a(l - p x )5 - Pi S u - <jS u - b( 1 - a)S u 


dS„ 


dt 


— Pl$v Pl$v 


e 1 s v 


dS,, 


dt 


®S U Pi $UC 5(1 ^Sy 




dt 


0 . 5(1 Pl)S v Pl^vnc 0 l^iSy 
— 5(1 0 )S u Pi^unc ~ vS U nc 


> (3.19) 


dE 

dt 


dt 

0iD ^vne T 5(1 c)S uc + vS unc (1 cci)E a ± E oc 2 E PiE 

dE T 


dt 


= a t E — PiE t — PiE t 


dl 


nT 


dt 


dI T 

-fa ~ 5(1 — <Xi)It ~ Pib ~ %ih 
— ^ 4(1 — 5)(1 — CC-P)!^ — mI nT — PiI n T — YI nT 


Let 


Then, 


Hi (t) — S(t) + S v (t) + S u (t) + S vc (t) + S uc (t) + Sync (0 + Sunc (0 

+ E(t) + E T (t) + I T (t) + I nT {t) (3.20) 

Wi(0 = 0i - PiS - a( 1 - p x )5 - PiS + 

+Ti[l ~ (5 + + S u + S vc + S uc + S vnc -1- S unc + E + E t + I T + I nT )] 

-\-piS PiS v PiSy < 3 : 5(1 Pi)S v + o(l Pi)5 Pi$u @S U 5(1 <3")5 U 

+ 01 Sy plSyc ^l^vc 4” <3"S W Pl^uc 5(1 <-)5 wc T 05(1 Pi)Sy PlSynC 01^1 Sync 

+5(1 <^)5 W Pl^unc vS unc + 0i3"i Sy n c + 5(1 c^S uc + vS unc (1 &i)E (%iE 

—a 2 E — p t E + a x E — PiE t — PiE t + 5(1 — a.i)I T — p 1 I T — ^ I T 
+ /c 4 (l — 5)(1 — CCi)I n r — m ^nT — Pl^nT ~ Y^nT 


Here it is important to note that in the absence of the disease, 

a 1 E, (1 — a 1 )E,a 2 E, 5(1 — a 1 )I T ,^ 1 l T ,Yl n T>^'iO- — 5)(1 — cti)I nT ,ml n T will be zero. 


Therefore, we have 


tfr(t) = 

01 (.PlS + PlSy + PlSu + PlSyc + PlSuc + PlSync + PlSunc + 
+ y 1 (l-H 1 ) 


PiE + PiE t + pJr + pJnr) 


Hi(t) = 0i - p 1 H 1 + y 4 ( 1 - // 4 ) 
ffi(t) = 0i + Ti - (Pi + Yi) h i 
Hi(t) + (pi +7 i)Hi < 0i +Yi (3.21) 


Now, we integrate both sides of (3.21) using integrating factor method (Kar and Jana, 2013; Birkhoff and 
Roffa, 1989). 
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That is, 

Hi +g 3 (t)dt = f(t) 

=^> H t < e -Js 3 (t)dt e la 3 mtf^ dt + c 2 ' 

Where g 3 (t) — + y 1 and f{t) — t/q + y 1 . The integrating factor is 

r(t) = e^ a3 ^ dt — e^ ai+Yl ^ dt = e^ 1+rY>t . 


Integrating (3.21) using the value of the integrating factor r(t) = e 


— pOi+ZiP 


we get 


tfi(t)<^y J f(t).r(t)dt + c 2 

Hl(t) ~ e (J+n)t [/ + rJ- e (fll+Vl)t dt + c 2 


< 


e (Pi+7i)t 

Htit) < 


</>i + Y i 


_ e (Mi+7i )tdt + c 7 


P + Ti 

^ 1 + Yl m e ^i+Yi)t dt + C2 

LMi + Y 1 

... H ( t ) < mlh + ce -(M+Yi)t 
Pi + Ki 


>-Oi+7i)t 


(3.22) 


Where c 2 is the constant of integration and if we let t -> qo, we have 

Hi(t) — ——— = [S + S v + S u + S vc + S uc + S vnc -I- S unc + E + E T + I T + I nT ] 
Pi + 7i 


(3.23) 


But 


5 = 


ViQti + OJfai+Pi + ke) 


(Pi + Pi + ^ 5 )(Pi + #i)(Pi + Pi + k 6 ) - 9 1 f 1 p 1 


We observe from the dynamic described by (3.1), (3.7) and (3.23) that the region 

f (s, S v , S u , S vc , S uc , S vnc , S unc , E, E t , I t , I nT ) G R+ 1 : 

^i(Pi + 0i)(Pi + /?i + K) 


Ti 


S < 


(Pi + Pi + k 5 )(pi + #i)(Pi + Pi + k 6 ) — Pi/?iPi 


[S + S v + S u + S vc -I- S uc + S vnc + S unc + E + E t + I T + I nT ] < 


V 1 +Y 1 
Pi + 7i^ 


is positively invariant. In other words, the system is well posed mathematically since it has a positive invariant. 
Therefore, for the initial starting point H 1 E M+ 1 , the trajectory lies in T^ Thus, we focus our attention only on 
the region T v 

3.9: BASIC REPRODUCTION NUMBER 

Definition 1': Basic reproduction number of an infectious disease model denoted by R 0 is the average number 
of secondary infections caused by a single infectious individual during their lifetime (Driessche and Watmough, 
2002). When the basic reproduction number is less than one, (P 0 < 1), the disease free equilibrium point is 
locally asymptotically stable. But if the basic reproduction number is greater than one, ( R 0 > 1), then the 
disease free equilibrium point is unstable. 

We shall calculate R 0 using the next generation matrix by starting with those equations of the model that 
describe the production of new infections and change in state among infected individuals. We refer to these set 
of equations as the infected sub system and linearize them about the disease free steady state. From equation 
(3.24) we derive the infection subsystem as transmission part and transition part respectively. 


@ IJTSRD I AvailableOnline@www.ijtsrd.coml Volume-2 I Issue - 5 IJul-Aug2018 


Page: 1663 


























International Journal of Trend in Scientific Research and Development (UTSRD) ISSN: 2456-6470 


(IS, 


dt 


r = PiS - Ml S v - PiS v - a 5 ( 1 - p t )s v 


dt 


= a( 1 - pjS - mA - oS u - b( 1 - a)S u 




dt 


PlSv Ml ‘-’vc 


dS, t 


dt 


El^UC ^(1 ^~)Sy 


dS„ 


dt 

dS-u 


( 25(1 Pi)S v Ml Sync 01^1^1; 
— b(l — a)S u — Mi S unc t )S U nc 


> (3.24) 


dE 

dt 


dt 

tpiT-i^vnc T d(l c)S uc + vS unc (1 a^E (XiE cc 2 E Mi E 

dE T 


dt 


— a t E — Mifi'r — Pi^t 


dl 


nT 


dt 


dI T 

~^T — ^(1 — a i)h ~ Pih ~ %ih 
= k 1 (1 - d)(l - a^)I nT - mI nT - Mi I nT ~ YKt 


Let M — [S,S V ,S U ,S VC ,S UC , S vnc , S unc , E, E r> I T , I nT ] r where T denotes transpose then we have the 
linearized infected sub system as transmission (F) and transition (V). 

M = (F + R)M (3.25) 

The matrix F corresponds to transmission part and it describes the the number of ways that new infections can 
arise. The matrix V corresponds to transition part which describes the number of ways that individuals can 
move between compartments and it includes removal by death or acquisition of immunity. Therefore, from 
(3.24) we obtain the matrices 


F := 


0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 (1 -a)b 0 0 0 

0 0 


0 (1 - c)5 <j> jT J 


0 

0 

0 

0 


0 

0 

0 


0 

0 

0 


0 

0 

0 

0 


0 

0 

0 


0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

v 0 0 


0 0 ttj 0 

0 0 0 


0 

0 

0 

0 

0 

0 

0 


0 0 0 ( 1 - a.) 
0 0 0 0 


0 

0 

0 

0 

0 

0 


0 

0 


(1 - d)|l - ajk. 
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0 + m+yj 


,-1 


o^-b^-ob-n tij '^Sbc-ocMj 




0 




i 


v -^ 
0 0 


0 

1 




-l 


0 0 


0 0 


0 0 


0 0 


0 0 


0 0 


0 0 


P-^ 




0 0 


y-m-u 
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FV - 1 • 


a < p.-vT 1 4, 


0 

0 

0 

0 

b-(c-p 

b+Mj 

S-g-(c-I) 


&(c-0 


• 


a.-u. +|Jli +a.T.4» 1 +n 1 -T.4 1 S-4, +b-u. -t-o-fc -r. 4 . ‘-o-b-c-o-c-4. 

0 0 0 

0 0 0 

0 0 0 


o-c-C+4, T,4,+4- v+4. 


a. 


0 0 


4 1 +a,+l 
0 0 0 


0 

a. -1 


- 4 - 


0 0 


0 

0 

0 

0 

0 


0 

0 

0 


(d-0ios-Uc 4 

Y+nw- 4 , 


By properties of determinant, we have that 

> If every element in a row (or column) of an n x n (square) matrix is zero, then the value of the determinant 
will be zero. 

> If the corresponding elements are equal in two rows (or column) of an n x n (square) matrix, then the 
determinant will be zero. 

Since the matrix above satisfies these two properties, we therefore state that the 

|FR _1 | = 0 

Here, we can now obtain the basic reproduction number R 0 from the trace and determinant of the matrix by 
letting FV -1 = G. 

Let 

R 0 — m(G) = - [trace (G) + Jtrace(G ) 2 — 4 det(G)^ (3.26) 

Observe that det(G) = 0 
So, 


Rn = 


V 


+ 


a 1 


v + Ml Ml + «2 + 1 


fc 4 (l — d)(l — ax) ^ 1 — a 1 


Mi + m + y 


Mi + ^1 


(3.27) 


Thus, if R 0 < 1, then the disease free equilibrium is stable and HBV infection can be controlled using the three 
control measure. On the other hand, if R 0 > 1, then the disease free equilibrium is unstable and the HBV 
infection increases (persists’) in the population. 


3.10: ENDEMIC EQUILIBRIM POINT 

This is the point where the disease cannot be totally eradicated but remains in the population. Let E **be the 
endemic equilibrium point. Therefore, for the disease to persist in the population, E** must not be equal to zero 
at equilibrium. To calculate the endemic equilibrium point, we have 

77** r r** r 1 * r** r** r** r** r** 77* 77* 1* j* -f- H 

C, , o v , o u , o vc , o uc , o vnc , o unc , CL , Cy, If, l n j, t\ J f U 
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If we solve equation (54) simultaneously, having in mind that E** 4 Owe have that; 

d( 1 - ct^Ij - Mr It ~ Si It = 0 
(d( 1 - «r) - (Mi + Si))I *t = 0 (3.28) 


k,{ 1 - d)( 1 - aM T - ml* nT - Mr I* T ~ Y^r = 0 
(/c 4 (l - d)( 1 - a 4 ) - (m + Mr + y))/nr = 0 (3.29) 


7* — 
*T ~ 


Si It + yl*nT + - Mi R* - YiR* = 0 

^i+Yi)R*-a 2 E*-YlnT 

Si 


I* — 
l nT — 


(M 1 + Yi)R*-oc 2 E*-SiIt 


(3.30) 

(3.31) 


If we substitute (3.30) into (3.28) we obtain 

- «i) - (Mi 
(Mi + Yi)R* ~ a zE 


L KJ \ U ) VV VJULCll 11 

(d(l - a 4 ) - (mi + ^i)) ^ 


(Hi + Yi)R*-a 2 E*- Y I^ 


Si 


= o 


j* _ 

l nT ~ 


Si 


(3.32) 


Substituting (3.31) into (3.29) we have 

1(1 - a 4 ) 

(Mi + Yi)R* - oczE* 


) iiilw yu j vv navu 

(^r(l - d)( 1 - a 4 ) - (m + Mi + y)) ^ 


(Mi + 7 i)R* - a 2 E* - 


y 


= 0 


r 

7? 


y 


(3.33) 


Observe that 
And 

Then, we obtain 


SiI*t + YInT + a 2 E* - Mi R* ~ YiR* = 0 

(d(l - aj ~ (Mi + Si))It = 0 
(/c 4 (l - d)(l - a 4 ) - (m + Mi + y))/nr = 0 


a 2 E* — Mi R* — y x R* — 0 
a 2 £* 

R* - —-- 

Mi + yi 

a x E* - Mi E* t - p x E* T = 0 


Et = 


a x E* 


(Mi + Pi) 


(3.34) 


(3.35) 


01 *iSvnc + 5(1 - c)S£ c + vSunc ~ (1 ~ Ui)E* - a x E* - a 2 E* - Mi E* = 0 
E , = (pl^lSync + 5(1 - c)S£ c + vSu nc 36 ^ 


Mi + a 2 + 1 

b(l-a)S* u - l i 1 S* unc -vS* unc = 0 


r** _ 

^unc 


5(1 ~ (t)S„ 
Mi + v 


(3.37) 
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(3.38) 


Let 


S* — 
°vnc 


a s(l Pi) $v Pl^vnc < fil T lSvnc ~ 0 
a 5 (l-0i)S v * 


Ml + 01 T 1 

aSZ ~ Mr S u * c ~ 5(1 - c)S u * c = 0 


* 

uc Mi + 5(1 - c) 


n* _ 

^VC 


PlSv 

Mi + ^1 


p 1 s;-ii 1 s* vc - 0^ = o 


a(l - p 1 )S* - Mi5* - aS* - 5(1 - <r)S* = 0 
a(l — Pi)S* 


S* — 

" ! Mi + o' + 5(1 — o') 


c* — 

‘-'P 


Pl 5* - Mi S v * - P 1 S; - a 5 ( 1 - /?,) 5* 
PiS* 


Mi + 0i + a 5(l _ Pi) 

xpi — Pl S* — a(l — Pl )5* — MiS* + 0i5^ c + Vi^t + 7i^* = 0 

PiS* 


(3.39) 

(3.40) 

(3.41) 

(3.42) 
(3.43) 


* - - a (1 - Pl )5- - ^ ^ ) 




+ 


V 1«1 


(Mi + Pi) 


0 lT 


05(1 — Pi) 


PiS* 


+5(1-c) 


1 V Mi + 0iTi / Vmi + Pi + a 5 (l - p 

a \ ( a(l — Pl )S 


+ 05(1 _ Pi)' 

a) ' 


_ \ ( o(l 

- c)/ \Mi + o 




+u 


+ 5(1 - cr) / 
a(l- Pl )5* \ 




+ 


r i«2 
Mi + 7i 




Mi + 5(1 
5(1-oK 

Pl + V J \pi + <7 + 5(1 — < 7 ), 
/ g 5 (l -0i) \ / Pi 5* _\ 

1 1 \ Mi + 0i T i / 'Mi + Pi + a 5(l _ 0i)' 

, xn _ , ( O \ ( _a(l - Pi)S* 

( C Hm! + 5(1 - c)J VM! + O 
^5(1 -ff)5* 


( _i_) 

\Pi + cc 2 + 1 / 


/ 


\ 


+u 


+ 5(1 - a), 

_ a(l ~ Pl )5* \ 

Pi + v )\pi +a + b(l-a)J / 


Mi + « 2 + 1 


= 0 


c = 


B — — 


01 T 1 


A = (fii + a( 1 - Pl ) + Mi) 
0i / Pi 


Pi + 61 \pi + 0i + a 5 (l - /? 

( 0^5(1 Pi ) 


a) 


(Mi + a 2 + 1) \ \ Mi + 0 iTl 7 Vmi + Pi + a 5 (l ~ 0: 


Pi 


y) 


D — — 


] = 


o 


a( 1 - Pi) 


5(1 ~ c) 

(Mi + 1) \ V P i + 5(1 - c)7 \ P i + o + 5(1 - a) 


'b(l-a)S* 


a(l- P i) 


(Mi + a 2 + 1) V \ Pi + v J V pi + o + 5(1 — a) 
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*/>i 


S* = 


A+B+ (n7^) ic+D+n+ (.KTr:) ic+D+n 


(3.44) 


3.11: SENSITIVITY ANALYSIS OF THE BASIC REPRODUCTION NUMBER OF HBV 
INFECTION 

The basic reproduction number R 0 is a function of eight parameters p 1; a 2 , v, a lt d, m, y and To control the 
disease, these parameter values must control R 01 such that it will be less than one(R 0 < 1). Therefore, as the 
parameter values are changed, we find the rate of change of R 0 . If we let g A — (p 1; a 2 , v, a lt d, m, y, ^ 1 ), then 
the rate of change of R 01 for a change in the value of parameter g 4 can be estimated from a normalized 
sensitivity index 

dR 0 R 0 

(3.45) 


zj° = 

9 4 


dg 4 9a 

Table 1: The Parameter values for Hepatitis B virus infection used to calculate R 0 . 


Parameters 

Description 

Value Range 

Reference 

P! 

Natural death rate 

0.0219 

Adu et al, 2014 

a 2 

Movement rate of non treated exposed (infected) 
individuals to recovered class 

0.88 

Abdulrahman et al, 2013 

V 

Movement rate of susceptible individuals who are 
neither vaccinated nor use condom 

0.002 

Assumed 

C*! 

Treatment rate of exposed (infected) individuals 

0.7 

Kamyad et al, 2014 


Movement rate of treated infectious individuals to 
recovered class 

0.98 

Kimbir et al, 2014 

d 

Treatment rate of infectious individuals 

0.98 

Kimbir et al, 2014 

m 

Rate of disease induced death 

0.0007 

Abdulrahman et a, 2013 


Movement rate of non treated infectious 

0.508 

Abdulrahman et al, 2013 

y 

individuals to recovered class 


To calculate the value of R 0 , we use the parameters as stated in table 1. 

v a x /c 4 (l — d)(l — a x ) 1 — a t 

0 v + pi Mi + a 2 + 1 Pi + m + y 9i + ^ 
0.002 0.7 0.006 0.3 

--+-+-+-= 0.76262 

0.0239 1.9019 0.5236 1.0019 


The normalized sensitivity index of the basic reproduction number with respect to p 4 , a 2 , v, a x , d, m, y and ( r 1 
is given by 

ryRg _ dR 0 R 0 _ Pi (R 0 ^ 

4 v 


dv ' v (v + p 4 ) 2 




14619.311 


v 


+ 


a i 


(v + Pi) 2 (pi + a 2 + l) 2 
r 0 _ dRp R 0 ( 1 


+ 


7“ 0 — —l _ — 

ttl da 1 ' a x 


7 r 0 _ dRp Rp 

111 ~ dg, ‘ P! 
/c 4 ( 1 - d)( 1 - 

(p 4 + m + y) 2 
1 


1 a 1 \ /Rq\ 

+ 07T^)(v) = - 38 - 9114493 


- -- - ——) (—) = -0.39687895 

Pi + a 2 + 1 Pi + <(i pi + m + y I Va 


rj R 0 _ dRp RP _ 

“ 2 da 2 ' a 2 


CCi 


(Rp 

(Pi + a 2 + l) 2 \a 2 


O- 


-0.31895975 
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yR 0 _ dRo Ro ^ 4(1 — a i ) (Rq n 
~~ ~d ' 


a dd 

y R 0 _ 

dm m 


&)- 


^! + m + 7 
/c 4 (l - d)(l - a x ) /R 0 


-0.44586536 


dRo Ro 
dy ' y 


^o\ _ 
m) 

(-) = 
(p-^+m + y ) 2 \y J 


(jj . 1 + m + y ) 2 
k 4 ( 1 - d)( 1 - a 4 ) (R q 


-23.8430637 


-0.03285462 


7 r q _ dRo Ro 


(1 ~ a 1 ) 

(Mi + £l ) 2 


= -0.32712716 


Table 2: The effect of the parameters on R 0 . 


Parameters 

Value Range 

Effect on R 01 

V 1 

0.0219 

-38.9114493 

a 2 

0.88 

-0.31895975 

V 

0.002 

14619.311 

a A 

0.7 

-0.39687895 

fi 

0.98 

-0.32712716 

d 

0.98 

-0.445865 

m 

0.0007 

-23.8430637 

V 

0.508 

-0.03285462 


The sensitivity index Z{y) is positive and this shows that the value of R 0 increases as the value of v increases. 
The remaining indices Z , Z,Z (a 2 )< Z (^ 1 ), Z (y), Z (d) and Z(m ) are negative, indicating that the value 
R 0 decreases as p 1( a v a 2 , y, d and m increases. Actually, the effectiveness of control may be measured by 
its effect on R 0 . Once the parameters i 2 1 ,a 1 ,a 2 ,^ 1 ,y,d and m will reduce R 0 < 1, then it is curative if the 
reduction can be maintained. 

Therefore, the endemicity of the will be reduced if there is increase in the value of the following parameters; the 
rate at which exposed (infected) individuals obtain treatment (a 4 ), natural death rate of the exposed and 
infectious individuals (pj), spontaneous recovery rate by exposed individuals (a 2 ), disease induced death rate 
of infectious individuals (m), rate at which treated infectious individuals move to recovered class ( 7 ), rate at 
which infectious individuals receive treatment (d) and the rate at which the received treatment is effective (^j). 
This implies that these parameters can help in reducing the rate of HBV infection over time and if it is 
maintained, the transmission of the disease may decrease, causing the cases in the population to go below an 
endemicity threshold. 


IV. Analysis of Result 

We consider the parameter set A= {xp, a, b, a s ,a, 5(1 — c), v, a v a 2 , d,k 4 , ,y,m,p 1 ,y 1 ,6 1 ,p 1 }; 

some of these parameters are taken from published articles and some are assumed with feasible values. More 
so, we consider V= {S 0 > S u o,S vc o,Suco. S V nco> Sunco> E 0 ,E TQ ,l T 0 ' l nT 0 > as the initial condition for 
the simulation of the model. We compare the results; having no control and having both vaccine, condom and 
treatment controls. With the parameter values in table 1, the system asymptotically approaches the endemic 
equilibrium E** — (8.62693005,0.86183812,0.09575979, 

0.95879491, 0.3836662,1.5315378,0.52086904,0.073104, 0.073104, 0, 0,0.12663685). 

The basic reproduction number is given as; 

_ v a-! fc 4 (l-d)(l-a 1 ) 1 -a x 

An —- 1 - 1 - 1 -— 

17 +/^ / i 1 +a 2 +l [d 1 +m+y Mi+si 

and the value of the basic reproduction number is R 0 — 0.76262 which is less than one. Thus, since R 0 < 1, 
then the disease free equilibrium is stable and HBV infection can be controlled in the population using the three 
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control measure. This confirms the result of the sensitivity analysis of R 0 in table 2, where d (the rate at which 
the infectious individuals are treated), a 1 (the rate of the exposed individuals who obtain treatment) and a 2 (the 
rate of exposed individuals who spontaneously recover from the infection) are all negative and it shows that 
when all these parameters increases, the average number of secondary infections caused by a single infectious 
individual during their lifetime (I.e, R 0 ), will decrease. This implies that with effective control measure, the 
future number of HBV infection cases will decrease as the parameters a 1 ,a 2 ,fj. 1 ,^ 1 m,d,y increases because 
even though individuals remain infectious for a longer period and consequently can infect more individuals, the 
rate of progression to disease decreases. 
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Figures (2a - 2c): Numerical simulation of the basic reproduction number R 0 using different rate 
of d, a t and a 2 , where and a 2 are respectively represented with b4 and b5. 


The effect of the three control measures on the disease dynamics may be well understood if we consider Figures 
2 a, 2 b and 2 c which explains how control reproduction ratio R 0 evolve with different rate of d, a 1 and a 2 where 
a 1 and a 2 are respectively represented with b4 and b5 in figure 2b and 2c above. It is seen that condom, 
vaccination and treatment reduces the value of R 0 effectively. But an integrated control works better than either 
of the control measures. 
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Figure (3ai - 3axii): Numerical simulation of the model system (3.1) when there is no control with time ranging 

from 0 - 6 months. 
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Figure (3bi - 3bxii): Numerical simulation of the model system (3.1) when there is no control with time ranging 

from 0 - lyear. 
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Figure (4ai - 4axii): Numerical simulation of the model system (3.1) when there is control with time ranging 

from 0 to 6months. 
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Figure (4bi - 4bxii): Numerical simulation of the model system (3.1) when there is control with time ranging 

from 0 to lyear. 


Here, we numerically solve the control system (3.1) 
which consists of twelve differential equations. We 
have plotted susceptible individuals, vaccinated 
individuals, non vaccinated individuals, vaccinated 
individuals who use condom, vaccinated individuals 


who does not use condom, non vaccinated individuals 
who use condom, non vaccinated individuals who 
does not use condom, Exposed (infected) individuals, 
Exposed but treated individuals, infectious but treated 
individuals, infectious but not treated individuals and 
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recovered individuals, by considering values of the 
parameters in table 5 and we simulated the system at 
different level of p lt /?j o, a 1 and d. 

From figures (3ai - 3axii) and (3bi - 3bxii), we see 
that the number of both the exposed and the infectious 
population increases when there is no control. Also 
from (4ai - 4axii) and (4bi - 4bxii), we observe that 
the application of control (vaccination, condom and 
treatment) gives a better result than when there is no 
control because the population of the exposed and the 
infectious individuals is reduced. This implies that 
effective vaccination with appropriate use of condom 
(condom efficacy) is able to reduce the density of the 
exposed (infected) individuals, while application of 
treatment is able to reduce the density of the 
infectious individuals. Therefore, since our main 
purpose is to reduce the density of the exposed and 
infectious individuals and increase the density of the 
removed with or without permanent immunity, we can 
conclude that the numerical simulation result shows 
that combination of vaccination, condom and 
treatment is the most effective way to control 
Hepatitis B virus infection. 

V. Summary and Conclusion 

To control HBV infection at least 96% Or more of 
newborns should be vaccinated and effort should be 
made to vaccinate more than 50% of the sexually 
active susceptible individuals who should also be 
encouraged to use condom during sex as these 
combination will go a long way in reducing the 
density of the exposed (infected) individuals. Also, 
treatment strategies are needed to control the spread 
of HBV infection in the population since control of 
HBV infection implies reducing the susceptible and 
exposed by vaccination and reducing the infectious by 
treatment. 
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